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ABSTRACT 

We investigate gas accretion onto a protoplanet, by considering the ther- 
mal effect of gas in three-dimensional hydrodynamical simulations, in which the 
wide region from a protoplanetary gas disk to a Jovian radius planet is resolved 
using the nested-grid method. We estimate the mass accretion rate and growth 
timescale of gas giant planets. The mass accretion rate increases with protoplanet 
mass for Mp < Mcri, while it becomes saturated or decreases for Mp > M^i, where 
Mcri = 0.036 Mjup(ap/1AU)'''^^, and Mjup and ap are the Jovian mass and the or- 
bital radius, respectively. This accretion rate is typically two orders of magnitude 
smaller than that in two-dimensional simulations. The growth timescale of a gas 
giant planet or the timescale of the gas accretion onto the protoplanet is about 
10^ yr, that is two orders of magnitude shorter than the growth timescale of the 
solid core. The thermal effects barely affect the mass accretion rate because the 
gravitational energy dominates the thermal energy around the protoplanet. The 
mass accretion rate obtained in our local simulations agrees quantitatively well 
with those obtained in global simulations with coarser spatial resolution. The 
mass accretion rate is mainly determined by the protoplanet mass and the prop- 
erty of the protoplanetary disk. We find that the mass accretion rate is correctly 
calculated when the Hill or Bondi radius is sufficiently resolved. Using the oli- 
garchic growth of protoplanets, we discuss the formation timescale of gas giant 
planets. 
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INTRODUCTION 



Currently, more than 400 exoplanets have been observed. Most such planets are consid- 
ered to be gas giant planets. Although giant planets are preferentially observed, observations 
imply that gas giant planets, such as Jupiter and Saturn in our solar system, can be born 
around stars. Thus, it is important to understand the formation process of the gas giant 
planet. In general, gas planets are formed in the protoplanetary disk (or the circumstellar 
disk) around the protostar. However, there is a problems about the growth (or the gas 
accretion rate) of the gas planet, which is related to the resulting mass of a gas giant planet. 



According to th e core accretion scenario (jPerri fc Cameron! Il974l : iMizuno et al.l Il978 



Hayashi et al. Ill985l ). the protoplanet has a less massive hydrostatic gas envelope when a 
protoplanet core mass is less than Mcore ^ 10 M® where M® is the Earth mass, while the pro- 



accretion phase, to become a gas giant planet when M 


> lOMffi (e.g.. 


Mizuno et al. 


1978 




Mizund 


1980 


Stevenson 


1982 


Bodenheimer & Pollack 


1986 


Pollack et al. 1 


1996; 


Ikoma et al. 


2ood 


2001; 


Hubickvj et al. 


20051). Since the eas eiant planet acquires almost all of its mass 



in the runaway accretion phase, the gas flow and mass accretion rate in this phase are 
important to determine the protoplanet evolution and its resulting mass. Since gas accre- 
tion onto the protoplanet may be closely related to the formation of the cir cumplanet ar y 
disk and acquisition process of angular momentum (IMachida et al.l l2008l : iMachidal l2009l ). 
three-diinension al calculations are required for investigating the runaway accretion phase. 
Machidal (|2009[ ) showed that the circumplanetary disk is formed in a compact region near 
the protoplanet. Thus, we may have to resolve the present size of the gas giant planet to 
estimate the gas accretion rate onto the protoplanet. As well, to properly handle the outer 
boundary of the protoplanetary system (protoplanet and circumplanetary disk), the region 
sufficiently far from the gravitational sphere of the protoplanet, that is, the Hill radius, 
should also be included, since gas flows into the protoplanetary system from outside the 
Hill sphere. Therefore, there is a need to incorporate vastly different spatial scales, ranging 
from the planet current radius to the Hill radius. For example, the Hill radius of Jupiter 
('"H.jup = 5.4 X 10^^ cm) is about 700 times larger than Jupiter's current radius (7.1 x 10^ cm). 



in global t 


iree- dimensional calculations (e.g., Kiev et al. 2001: 


D'Aneelo et al. 2002. 2003: 


Bate et al. 


2OO3I Klahr & Kiev 20061: Dobbs-Dixon et al. 200'' 


^ PaardekooDer & Mellema 


20081: Fouchet & Mavei 20081: D'Aneelo & Lubow 2008h. However, the fine structures in the 



proximity of the protoplanet cannot be resolved in such calculations. A compact disk is 
formed in the range of r < 10 — 50 rp (IMachida et al.ll2008l : iMachidal l2009l ). Thus, the mass 
accretion rate may have to be derived by resolving the spatial resolution with at least ~ 10 r 



p- 
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It should be noted that the mass accretion rate obtained from a global simulation can be 
correct, if the mass accretion rate is determined solely by the global structure around the Hill 
sphere and the circumplanetary disk contributes little to gas accretion. Even in such a case, 
the mass accretion rate should be investigated in calculations with higher spatial resolution 
to confirm the validity of the accretion rate derived in simulations with coarser spatial reso- 
lution. There are only a few studies that report the mass accre tion rate onto the protoplanet 



with sufficiently high -spatial resolution in a local simulation ( iTanigawa &: Watanabd l2002l : 
Avliffe fc Bat^ boooh . However, in these studies, the gas flow from the region outside the 



Hill radius to the proximity of the protoplanet was not sufficiently discussed. Furthermore, 
the acquisition process of the angular momentum and formation of the circumplanetary disk 
are not discussed in these studies. 

In this paper, we focus on the mass accretion onto the protoplanet in runaway gas ac- 
cretion phase, in which the gas continues to collapse onto the protoplanet without additional 
heating by collision of planetesimals. After the solid core (or planetary core) formation, the 
formation process of the gas giant planet can be divided into two phases. The quasi -static 



envelope slowly contracts with a timescale of ~ 10 — 10 yr (e.g.. lIkoma et al. Il2000l ) when 
the protoplanet is less massive than Mcore ^ 1 ~ 10 M^, while the gas rapidly collapses 
onto the protoplanet when Mcore ^ 1 ~ 10 M^. We investigate the evolution of the proto- 
planetary system only in the runaway gas accretion phase, based on the results of a local 
simulation with a sufficiently high-spatial resolution using the nested-grid method, in which 
the region of ~ 5 — 20 th from the protoplanet is resolved with cells having the size of the 
radius of the present -day Jupiter. Since the acquisition pr ocess of the aii gular momentum 



jMachida et al.lboosh and circumplanetary disk formation (jMachida] l2009l ) around the pro 



toplanet were already investigated using this method, this paper will focus on gas accretion 
onto the protoplanet and gas flow in the proximity of the protoplanet. As a result of calcu- 
lation, we found that the mass accretion rate is correctly calculated when the Hill or Bondi 
radius is sufficiently resolved, and it barely depends on the thermal effect around proto-gas 
giant planet. The structure of the paper is as follows. §2 gives the model frameworks, while 
§3 describes the numerical methods used. The numerical results are presented in §4 and 
compared with the results of previous studies in §5. §6 discusses protoplanetary growth. 
The conclusions of this paper are presented in §7. 
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MODEL 



2.1. Basic Equations 



A local region around a pr otoplanet is considered using the shearing sheet model (e.g., 
Goldreich fc Lynden-Belll Il965[ ). in which the self-gravity of the protoplanet ary disk is ig- 
nored. In addition, no physical viscosity is included, and the numerical viscosity can be 
ignored because it is sufficiently small. Thus, an inviscid gas disk model is adopted. The 
orbit of the protoplanet is assumed to be circular in the equatorial plane of the protoplan- 
etary disk. Local rotating Cartesian coordinates with the origin at the protoplanet are set 
up, in which the x-, y-, and 2-axis are the radial, azimuthal, and vertical directions of the 
disk. The equations of hydrodynamics without self-gravity are solved [see, eqs. (l)-(6) of 



Machida et al. 2008 



For the gas, a barotropic equation of state is adopted (for details, see lMachidall2009l). In 



a loc al region, the protoplanetary disk has an almost constant temperature flHayashi et al 
19851 ). while the gas around the pr otoplanet, that is, the gas envelope, has a higher tempera- 



iysolj , wnne tne gas around tne pr otoplanet, tnat is, tne gas envelop 
ture than the protoplanetary disk (Mizunq et al. 1978[ lMizunolll980l: 



1986; Pollack et al. 1996; Ikoma et al. 



20001 



Mizuno et al. 



f ll978h " 



Bodenheimer fc Pollack 



studied the structure 

and stability of the envelope around the protoplanet, on the assumption that the envelope 
is spherically symmetric and in hydrostatic equilibrium. They also investigated the thermal 
evolution of the envelope, parameterizing the dust grain opacity, and d etermined the bound- 
ary between the isothermal and adiabatic regions. Using Figure 2 of iMizuno et al.l ( 1l978l ). 
the thermal evolution around the protoplanet is modeled as a function of the gas density, 
that is, using the barotropic equation of state, as 



Cs,0P 



tanh ( 

,Pcri 



+ Kp-^ tanh ( — 

. Peri 



where Cg is the sound speed, 7 is the adiabatic index (7 = 1.4), and the adiabatic constant 
K is defined as K = c^qP^^'^, where pcri is the critical density, wherein the gas behaves 
isothermally for p < pcri, and adiabatically for p > pcri- In this study, pcri = 00 (isothermal 
model), po, lOpo, 10^ po, and 10^ po (adiabatic model) are used, where po is the initial density 
on the equatorial plane. The hyperbolic tangent (tanh) function is used to smoothly connect 
the first (isothermal) and second (adiabatic) terms in Equation ([1]). The thermal evolution 
for different pcri is plotted against the gas density in Figure [T], in which the gas temperature 
is constant in the isothermal model (pcri = 00), while it increases gradually from the initial 
value at p ~ pcri in the adiabatic models (pcri = Po, lOpo, 10^ po, and 10'^ po). 



In the standard disk model (IHayashi et al. Ill985l ). at a Jovian orbit, the density is 
Po = 1-5 X lO^^^gcm"^, and the temperature is Tq = 123 K. Thus, for example, in a model 
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with Peri = 10 po; the gas behaves isothermally when p ^ 1.5 x 10~^°gcm~^, while it 
behaves adiabatical ly when p ^ 1.5 x 10~^°gcm~^. Comparing Figure [1] with Figure 2 of 
Mizuno et aLl (jl978[ ). the thermal evolution of the model with p^n = 10 po (Fig. □ broken line) 



corre sponds to that for a gas envelope around a proto- Jovian planet (Fig. 2 of iMizuno et al. 
19781 ) with a dust opacity of Kg = l.Ox 10 ~^ cm^ g~^. In t his rn o del setting, the critic al density 



Peri corresponds to the dust opacity Kg in lMizuno et al.l ( 1l978l ). iMizuno et al.l ( 1l978l ) adopted 
Kg = 1.0 X 10~^cm^g~^ as the most reliable parameter of a proto- Jovian planet, indicating 
that a more realistic gas temperature of the envelope is lower than that in the model with 
= 10 Po (the dotted line of Figure [1]). In our settings, because Kg = 1.0 x lO""^' 



cw? g 



almost corresponds to pcri = lO^po, we call the models having p^ri = lO^Po 'most realistic 
models'. It should be noted that the critical density pcri increases as the dust opacity Kg 
decreases. Thus, in models with pcri = lOpo, the thermal energy around the protoplanet may 
be overestimated. On the other hand, when the isothermal equation of state is adopted, 
the thermal energy around the protoplanet is obviously underestimated. Therefore, it is 
expected that the actual thermal evolution is located between models with pcri = po and 



Pc 



oo. 



2.2. Protoplanetary Disk Model and Scaling 



The ini tial settin g s are similar to lMiyoshi et al. I (119991 ) , iMachida et al. I (j2006bl ) , iMachida et al 



(2008 


), and 


Machida 


(2009) 



vo = (0, -3/21]pX, 0), 
where flp is the Keplerian angular velocity of the protoplanet 



GMr 



1/2 



(2) 



(3) 



where G is the gravitational constant, Mc is the mass of the central star, and ap is the orbital 
radius of the protoplanet. 

For hydrostatic equilibrium, the density is given by po = ctq exp (— z^/2/i^) / \^27ih, where 
ctq (= Po dz) is the surface density of the unperturbed disk. The scale height h is related 
to the sound speed Cghj h = Cs/flp. 



Using the standard solar nebular model (jHayashilll98ll : iHayashi et al. Ill985l ). the tem- 
perature T, sound speed Cg, and gas density pc,o can be described as 



T = 280 



Lcr 







1/4 



Qp y 

lAU/ 



-1/2 



K, 



(4) 
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where L and Lq are the protostellar and solar luminosities 



where /i = 2.34 is the mean molecular weight of the gas composed mainly of H2 and He, and 



1.9 X 10' 



cms 



(5) 



/a \ -11/4 
.1. 



po = 1.4 X 10"^ (^^j gem 



When Mc = 1 Mq and L = 1 Lq are adopted, the scale height h can be described as 

/a \ 5/4 
/i = 5.0 X lO^M— ^ cm. 
VlAU/ 

The inverse of the angular velocity Qp is described as 



(6) 



(7) 



VlAU 



yr. 



The basic equations can be normalized using unit time, f^p^, and unit length, h. The 
density is also scalable and is normalized using ao/h. Hereafter, the normalized quantities 
are expressed with a tilde on top, for examp l e, x = x/h, p = po/{ao/h), and i = tflp. 
Further details can be fou nd in iMachida et al.l (120081 ). A dimensionless description is given 
by Equations (14)-(19) of iMachidal (120091 ). The dimensionless quantities are converted into 
dimensional quantities using Equations (H])-®. The gas flow is characterized by two param- 
eters, the dimensionless Hill radius fn = rn/h, and the critical density pcri- The dimensional 
Hill radius is defined by 



Mp 
3iVt 



1/3 



(9) 



In this study, a Hill radius fn ranging from 0.29 to 1.36 is used. As a function of the orbital 
radius and the mass of the central star, the parameter fn is related to the actual protoplanet 
mass in units of Jovian mass Mjup by 



M 



0.12 



Jup 



IMo 



-1/2 



lAU/ 



3/4 



^3 



(10) 



For example, in the model with fn = 1.0, Op = 5.2 AU and Mc = 1 Mq, the protoplanet mass 
is Mp = 0.4Mj. In the parameter range of fn = 0.29 — 1.36, at a Jovian orbit (op = 5.2 AU), 
protoplanets have masses of O.OlMjup — IMj^p. Table 1 gives the dimensionless Hill (fn) and 
Bondi (fe) radii, the masses of protoplanets at Jovian (5.2 AU) orbit, the critical densities 
and sink radii for all models. Model names consist of two parts: the protoplanet mass at 
the Jovian orbit and the critical density. For example, model M001A3 has parameter values 
given as Mp = O.OlMj^p, p^ri = 10^ Po- 
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3. NUMERICAL METHODS 

3.1. Numerical Procedures 

The purpose of this study is to investigate gas accretion onto the protoplanet in three- 
dimensional simulations. However, given current CPU limitations, it is impossible to calcu- 
late the complete evolution of the gas giant planet with sufficiently high-spatial resolution, 
that is, the evolution of the planet from a solid core (~ 10 M^) with a thin gas envelope to a 
protoplanet that acquires a massive atmosphere (~ IMjup) cannot be performed. Thus, the 
mass accretion rate is derived using the following procedure: (i) With a fixed protoplanet 
mass, the evolution of the protoplanetary system until the gas fiow reaches the steady state 
(~ 100 orbits) is calculated, (ii) Then a sink is introduced, and the mass accretion rate that is 
derived is averaged over a further ~ 20 orbits, and (iii) Finally, steps (i) and (ii) are repeated 
by changing the protoplanet mass to give the growth rate of the protoplanet. 



3.2. Nested-Grid Method 



To investigate the formation of a circumplanetary disk in a protoplanetary disk, it is 
necessary t o cover a large dynamic spa tial scale range. Using the nested-grid method (for 
details, see lMachida et al. II2005I . l2006a] ). the regions near and remote from the protoplanet 
are covered with adequate resolution. Each level of a rectangular grid has the same number 
of cells (= 64 X 128 x 16), but cell width As{l) depends on the grid level /. The cell width is 
reduced by 1/2 with increasing grid level (/ — )■ / + 1). Eight grid levels (/max = 8) are used. 
The box size of the coarsest grid, / = 1, is {L^, Ly, L^) = (12,24,3), and that of the finest 
grid, / = 8, is {L^,Ly,L,) = (0.09375,0.1875,0.0234). The cell width in the coarsest grid, 
/ = 1, is As = 0.1875, and it decreases with As = 0.1875/2'"^ as the grid level / increases. 
Thus, the finest grid has As(8) ~ 1.46 x 10~^. The fixed boundary condition in the x- and 
5-direction, and the periodic boundary condition in the y-direction are used. 



3.3. Sink Cell and Smoothing Length 

In the finest grid (/max = 8), the cell width is As = 1.46 x 10"'^. In real units, when 
the protoplanet is located at 5.2 AU, the cell width corresponds to As = 5.7 x 10^ cm, or 0.8 
times the Jovian radius. The evolution of the protoplanetary system is calculated using a 
sink. In the fiducial models, the radius of the sink is Vsink = 3.53 x 10^^ or twice the Jovian 
radius at Jovian orbit (for details, see §3.41) . During the calculation, the gas from the region 
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inside the sink radius is removed in each time step. 

The smoothing length for the gravitational potential of the protoplanet is not explicitly 
used. For numerical calculations, the physical quantities are defined at the cell centre, while 
the origin (protoplanet 's position) is defined as the cell boundary. Thus, the region inside 
r < = -\/3As(/max)/2 has a uniform gravitational potential. At a Jovian orbit, is 0.7 
times the Jovian radius (fg = 1.26 x 10"^ or = 4.9 x 10^ cm). 



3.4. Convergence Test for the Accretion Rate 



To check the convergence of the mass accretion rate onto the protoplanet, the evolution 
of the protoplanetary system for different sink radii was calculated. The convergence of other 
quantities that change with the cell widt h or grid level w ith and without a sink were already 
investigated in lMachida et al.l (120081 ) and lMachidal ( l2009l ). Figure H] shows the mass accretion 
rate as a function of the sink radius for models M02ISS, M02ISM, M02I, and M02ISL. The 
mass accretion rate is derived based on the procedure outlined in §3.11 As listed in Table [H 
these models have the same Hill radius fn (or the same protoplanetary mass) and different 
sink radii rsink- Models have Mp = 0.2Mjup at Jovian orbit (op = 5.2 AU), and the isothermal 
equation of state is used. Model M02ISL has the smallest sink radius of fsink = 1.6 x 10"^ 
which corresponds to 6.3 x lO^cm (0.87 rjup), while model M02ISL has the largest sink radius 
of J^sink = 1-1 X 10"^ which corresponds to 4.3 x lO^^cm (6rjup). The protoplanetary system 
was calculated for ~ 100 orbi ts, for which it was confirmed that the gas flow achieves a 



steady state in t -C 100 orbits. iTanigawa fc Watanabe f|2002h a l so sho wed that steady state 



is reached after 1 — 10 orbits for a local calculation. iMachidal (120091 ) showed that the gas 
envelope and angular momentum around the protoplanet reach steady state after 1-10 orbits. 

Figure M shows that the models M02ISS, M02ISM, and M02I have almost the same 
accretion rate of dMp/dt ~ 0.18, for models, which have sink radii of fgink < 3.5 x 10~^ 
(1.4 X lO^^cm; 1.92 rjup). On the other hand, model M02ISL has a sink radius of fsink = 
1.1 X 10~^ (4.3 X 10^°cm; 6rjup), and has a slightly larger accretion rate of dMp/di ~ 0.35 
than the others. Thus, the difference in the accretion rate between M02ISS (0.2) and M02ISL 
(0.35) is only within a factor of 2. As a result, the accretion rate can safely be estimated 
when the sink radius is less than fsink < 3.5 x 10"'^. In the following, we show results adopting 
3.5 X 10-3 (1.4 X 10^°cm; 1.9 rj. 



'"sink 



J up; 



Since the gas inside the sink radius is removed at each time step as mentioned in §3.31 
the structure, that is, the protoplanet's atmosphere and envelope, inside the sink radius 
cannot be investigated. However, it is expected that the gas falling into the sink cannot 
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escape from the protoplanet and, hence, does not influence the results obtained. In the 
region where r < rsjnk, since the gra vitational energy greatly dominates the thermal energy 
( iMizuno et al.lll978l : iMachidal 120091 ) . the gas that falls into the sink cannot be pushed out 
by the gas pressure. In addition, the gas falling into the sink does not acquire additional 
angular momentum inside the sink radius (or near the protoplanet), because the angular 
momentum of the sys tem is only acquired by the shearing motion of the protoplanet ary disk 
( jMachida et al.ll2008l ). However, it is unknown how the gas trapped by the gravitational 
potential of the protoplanet reaches the surface of the protoplanet because the centrifugal 
force is expected to prevent the gas from further falling inside the sink. To understand 
the surface and structure of the protoplanet, the region inside the protoplanet, that is, 
inside r < rsink needs to be resolved. Such simulation would require an enormous amount 
of CPU time. In the calculations, the sink radius is rsink = l-9rjup that is smaller than 
the Roche limit = 2.5 Tp. Thus, it is expected that the gas inside r < rgink loses its 
angular momentum by tidal interaction and finally falls into the gas planet. In addition. 



Klahr fc Kleyl ( l2006l ) suggested that the radius of a young planet is about twice the current 



size of Jupiter. Therefore, when the sink radius is comparable to the size of the present 
planet, it is possible to safely estimate the gas accretion rate. 



4. Gas Flow and Mass Accretion Rate 



4.1. Spiral Pattern and a Circumplanetary Disk 

Figure E] shows the structure around the Hill sphere for each model at ~ 100 orbits 
(t ~ 640). The systems in Figure [3] are in steady state, because th e gas flow achieves a 
steady state in ~ 10 orbits (t ^ 6.3 — 63) as show n in §3.41 (see also iMiyoshi et al. I Il999 



Tanigawa fc Watanabd 120021 : iMachida et al.ll2008l ). Models in the figure have a parameter 
of Peri = lO^po; for which the gas temperature increases adiabatically for p > lO^po- The 
figures show that the glo bal structure in the models is almost the same as for the isothermal 
models (for details, see iMachidal l2009l ). In calculations using the isothermal equation of 
state, the spiral p atterns that are dis tributed from the upper-left to the lower-right region 
are observed (e.g.. lLubow et al. Ill999l ). The same patterns are seen also in Figure [31 in which 



the barotropic equation of state is used. This can be explained by noting that for regions far 
from the protoplanet, the gas density is not high enough, and the gas behaves isothermally. 
As well. Figure |3] shows that as the protoplanet mass increases, stronger shock waves appear 
and the spiral patterns become clearer. Furthermore, the density contrast between spiral 
and gap becomes stronger as the protoplanet mass increases. Thus, it can be concluded that 
the features observed in adiabatic models correspond well with those features observed in 
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isothermal models. 

To focus on the region near the protoplanet, the structures inside the Hill sphere are 
shown in Figure H] in three dimensions. The grid level of Figure S] is / = 6, while that of 
Figure [3] is / = 2. Thus, Figure H] is a 16 times enlargement of the central part of Figure [31 
In Figure HJ the region of p > lO'^po is plotted using the red constant density surfaces, while 
that for p > lO^po is plotted using the orange constant density surfaces. It should be noted 
that the orange and red surfaces do not appear in model M001A2, because this model has 
a very small part that has p > lO^po. The density distributions on x = 0, |/ = 0, and 2 = 
plane are projected onto each wall surface. Figure H] shows that the red and orange surfaces 
increase with the protoplanet mass, indicating that a massive protoplanet has a denser, or 
more massive, envelope. The orange surface flattens as the protoplanet mass increases, while 
the red surface in all mode ls except f or mo del M001A2 has a sufficiently flattened structure. 



Machida et al.l (120081 ) and iMachidal (120091 ) showed that the specific angular momentum of 
the envelope increases in proportional to j oc Mp when the protoplanet mass is smaller 
than M < Mj^p. Thus, the angular momentum also increases with the protoplanet mass. 
Therefore, the centrifugal radius increases with the protoplanet mass. Since the centrifugal 
force gradually affects the more distant region from the protoplanet as the protoplanet mass 
increases, the low-density region gradually flattens out owing to the centrifugal barrier. On 
the other hand, the red surfaces always have flattened structures, because the red surface (or 
higher density envelope) is located near the protoplanet that is inside the centrifugal radius 
even when the protoplanet is less massive Mp ~ O.OSMj^p. 



4.2. Mass Accretion Rate onto a Protoplanet 

Figure [5] shows the mass accretion rates onto the protoplanet for all models against 
the cube of the Hill radius, f^. With a fixed protoplanet's orbit, can be connected to 
the protoplanet mass using Equation ( fTOj) . As a reference, the mass accretion rate and the 
protoplanet mass at a Jovian orbit are given in the upper and right axes in Figure O The 
mass accretion rate and the protoplanet mass at the Jovian orbit can be converted into those 
at any orbit using the value of Op in parenthesis in the upper and right axes. In this paper, 
the mass accretion rate at the Jovian orbit will be used for convenience. 

As shown in Figure the mass accretion rates increases rapidly for protoplanetary 
mass in the range of Mp < 0.2Mjup. The accretion rates for isothermal and adiabatic 
models are almost the same in this range. The accretion rate for each model has a peak 
at Mp ^ 0.2Mjup, then it gradually decreases in the range of Mp > 0.2Mjup. As well, in 
this range, the accretion rates depend slightly on the equation of state used. The difference 
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increases as the protoplanet mass increases. For Mp > 0.2Mjup, the decrease of the mass 
accretion rate in models with a harder equation of state, that is smaller pcri? is less than 
in models with a softer equation of state, that is, larger pcri- For example, the barotropic 
model of Pcri = 10 po has accretion rates of Mp = 4.7 x 10~^Mjupyr~^ at M = 0.2Mjup, and 
4.1 X 10~^Mjupyr~^ at IMjup. On the other hand, the isothermal model shows a steeper 
decrease of the accretion rate than the other models. The isothermal model had accretion 
rates of Mp = 4.3 x IQ-^Mj^p yr-^ at Mp = 0.2Mjup, and 2.3 x 10~^Mjup yr"^ at M = IMj^p. 
Thus, the difference between isothermal and barotropic models is not so large even when the 
protoplanets become massive. When the protoplanet has a Jovian mass IMjup, the accretion 
rate in model with p^ri = Po is only 2.4 times larger than that for the isothermal model. 

As shown in §2.11 it is expected that the actual mass accretion rate is located between 
the isothermal and the adiabatic model with pcri = Po- Around the protoplanet, the thermal 
energy is underestimated in the isothermal model, while it is overestimated in the model 
with Peri = Po- We fitted the mass accretion rate by the solid line in Figure [5l in which, 
for convenience, it is fitted as a constant in the range of Mp > 0.125Mjup. It should be 
noted that, although the accretion rate in each model gradually decreases in this range, the 
rate of decrease is sufficiently small especially in adiabatic models. As well, it can be noted 
that, since the protoplanetary evolution in the range of Mp > IMj^p was not calculated, the 
accretion rate may not apply in the range of Mp > IMjup. Using dimensionless quantities, 
the mass accretion rate is described as 



dMr. { 0.83 ff?, )3/2 for fl < 0.3 



dt 1 0.14 for f3 > 0.3. 

In Equation (fTTI) . the mass accretion rate for t\ > 0.3 {Mp = 0.14) almost corresponds to 
that of the average (Mp = 0.16) in most realistic models having p = lO^po for fn > 0.3. 
Since the protoplanetary evolution was calculated using a local simulation, it is possible to 
convert the mass accretion rate from the fixed orbit to any arbitrary orbit, arbitrary density, 
and temperature of the protoplanetary disk model. Using Equations (jl])-®, the dimensional 
mass accretion rate can be described as 



dMp 
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When the central protostar has a mass of M 

standard model (i.e., MMSN disk model; iHayashi et al. Ill985[ ) with surface density Sq is 



adopted in Equation (fT2l) . the accretion rate at any orbit can be described as 



dMr, 
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(13) 



where S is arbitrary surface density at any orbit. Equation (1131) indicates that the growth 
timescale of the Jovian mass protoplanet at a Jovian orbit is ~ 10^ yr. We will present a 
detailed discussion of the mass accretion rate and the growth rate of the protoplanet in §5. 
The following sections show the gas flows at large ( §4.31) and small scales ( §4.41) around the 
protoplanet to confirm the relationship between the accretion rate and the circumplanetary 
disk. 



4.3. Gas inside the Hill Sphere 



This section examines the large-scale gas flows in order to determine gas accretion from 
outside the Hill sphere. Each top panel in Figure [6] plots the gas streamlines falling into 
the protopl anet. Figure shows only models w ith pcri = lO^po, which is the most realistic 
parameter (IMizuno et al.lll978l : lMachidall2009l ). There are no differences in the streamlines 
between models with different equations of state, because the gas outside the Hill sphere, 
that is, at a large-scale, has a low density and behaves almost isothermally. In the figure, the 
streamlines are inversely integrated from the region inside the Hill sphere. Thus, only the 
streamlines accreting onto the protoplanetary system are plotted. Colors of the streamlines 
indicate the velocity in the vertical direction (w^). For example, the gas rises upward in blue 
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and green parts of the streamline, while it goes down in yellow and red parts. The Bondi, 
which is shown using a black circle, 



(14) 



the doublewide Bondi, shown with a black-dotted circle, the Hill, shown with a white circle, 
and the doublewide Hill, shown with a white-dotted circle radius, are also plotted in the 
figure. Each top panel in Figure [6] shows that the gas flow falling into the protoplanetary 
system is controlled either by the Bondi or the Hill radius. When the Bondi radius is smaller 
than the Hill radius (models M001A2 and M005A2), the gas flows into the protoplanet only 
in the narrow band of tb ^ |x| < 3 tb- The bandwidth of streamlines in the region far from 
the Hill sphere almost corresponds to the Bondi radius (~ tb, see the lower left corner in 
the top panel of Fig. [6]) . When the protoplanet mass increases and the Bondi radius exceeds 
the Hill radius (M04A2, M08A2, and M1A2), the gas flows into the protoplanet in the band 
of '"H ^ ^ 3rH, in which the streamlines in the region far from the Hill sphere have a 
bandwidth of the Hill radius (~ rn). It should be noted that in Figure El the grid-scale 
is different in each of the panels. Thus, it can be concluded that the gas flows into the 
protoplanetary system occur only in the region given as 



< a; < 3 r. 



race = mm(rB,rH 



(15) 



It was also shown by iD'Angelo fc Lubowi (120081 ) that the mass accretion rate is controlled 
by either the Bondi or the Hill radius. 

The colors of the streamlines imply that the gas rises slightly in regions far from the 
protoplanet (see, for example. Models M04A2, M08A2, M1A2), then it gradually goes down 
as the streamlines approach the Bondi or Hill sphere. Near the Bondi or Hill sphere, the gas 
rises strongly upward near the shock fro nt and then rapidly fa lls into the Bondi or Hill spheres 
in the vertical direction (for details, see lMachida et al.ll2008f). The s i milar f eatures of the gas 



flow ( or str eamlines) around the Hill sphere are seen in iKley et al. 
(I2OO3I ) and lPaardekooper fc Mellemal (|2008l ). 



fl200lh . ID'Angelo et al. 



For each bottom panel in Figure [6l the plane is plotted, in which the color indicates 
the mass flux p\v\ per unit area. In the plane, only gas passing through the gray region (or 
gray points) can reach the protoplanetary system (compare the gray regions in each bottom 
panel with the streamlines in each top panel). Most of the gray points are distributed 
in the region far from the equatorial plane, indicating that the gas present in the midair 
flows preferentially into the protoplanetary system, while the gas near the equatorial plane 
barely flows into it. Although the gas near the equatorial pla ne also flows into the Hill 
sphere, it flows out from the Hill sphere afterwards, as shown in iMachida et al.l fl2008[ ) and 
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Paardekooper fc Mellemal (120081 ) . Figure [7] shows the streamhnes for model M04A2, in which 
the streamhnes are integrated from the area of fn < x < 2fH and < z < 1.5 h on the 
y = —6h plane, that is, the region denoted by the black squares. It should be noted that the 
streamlines in Figure E] are mverse/?/ integrated from the region inside the Hill sphere, while, 
in Figure [TJ they are integrated along the flow from inside the black square. The color of each 
streamline indicates the mass flux p\v\ per unit volume at each mesh point. The gas flow 
through the black square in Figure [7] tracks three different paths: (i) the flow is attracted 
toward the protoplanet and returns to a similar orbit as the original orbit (the pass-by 
region), (ii) the flow is trapped by the protoplanet and falls into the protoplanetary system 
(the planetary atmosphere region), ar id (iii) the flow turns round inside the Hill sphere and 



goes back out (the horseshoe region). iLubow et al 



in a certain band (see. Figure 4 of iLubow et al. 



(119991) showed that the flow distributed 



19991) accretes onto the p rotoplanetary 



system in a two-dimensional simulation (see also. iTanigawa fc Watanabell2002l ). However, in 
three-dimensional simulations, only a part of the flow in the band (or bundle) i n the radial 



direction fl ows into the protoplanetary sys tem, as shown in Figures [6] and [71 iKley et al. 



(I2OOII ) and lPaardekooper fc Mellemal (120081 ) pointed out that the mass accretion rate derived 
in three-dimensional calculations is much smaller than that derived from two-dimensional 
calculations. 

As shown in Figure |6l the bandwidth that limits the gas flowing into the protoplanetary 
system is proportional to mm(rs, th). The Bondi radius is proportional to oc (eq. [T4|) . 

1 /3 

while the Hill radius is proportional to oc Mp (eq. When the Bondi radius is smaller 
than the Hill radius (re < th), the bandwidth (~ re) and accretion rate rapidly increase with 
the protoplanet mass because the mass dependence of the Bondi radius is strong (re oc M^). 
On the other hand, when rn > tb, the bandwidth (~ th) and accretion rate maintain 
an almost constant value, because the mass dependence of the Hill radius is considerably 

1 /3 

weaker (rn oc Mp' ). It should be noted that the Bondi radius equals the Hill radius when 
Mp = 0.023 (ap/lAU)^/^Mjup. Thus, as shown in Figure O it is possible to quahtatively 
understand the protoplanet mass dependence on the accretion rate: the accretion rate rapidly 
increases when re < th, while it maintains an almost constant value when re > th. However, 
the power of the accretion rate (Mp oc M^'^ for tr < ^h, see Fig. derived in our calculations 
is slightly smaller than the Bondi solution (IBondil 119521 ): 



ciMp _ AttG^M^po 



dt 



oc M' 



(16) 



--5,0 



This discrepancy can be attributed to the fact that spherically symmetric flow in an isolated 
system is assumed in Bondi accretion. However, in the protoplanetary system, the gravita- 
tional sphere is limited to the Hill sphere, and the angular momentum that was acquired 
from the shearing motion in the protoplanetary disk can affect the gas accretion. Thus, it is 
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natural that the power of the accretion rate in the protoplanetary disk does not completely 
match with that obtained from the Bondi solution. In fact, the accretion rate derived in 
this simulation is two or three orders of magnitude smaller than the Bondi accretion rate. 
As has been mentioned, the Bondi radius is useful in representing the gas flow when the 
protoplanet is less massive, indicating that the thermal pressure cannot be ignored when 
the thermal energy dominates the gravitational energy inside the Hill sphere. Note that, in 
the isothermal model, the thermal energy d ominates the gravitat ional energy inside the Hill 
sphere when Mp < O.OSMjup (for details see lMachida et allboosl ). 



4.4. Gas Flow around a Protoplanet 

In this subsection, we consider the gas flow patterns near the protoplanet. Figure [8] 
shows the density distributions and streamlines in the neighborhood around the protoplanet 
for models M005A2 (left) and M1A2 (right). The streamlines are inversely integrated from 
the high-density region (yellow constant density surface). The left panel in Figure [H] shows 
that the gas spirals into the protoplanet in the vertical direction. On the other hand, in the 
right panel, a part of the gas accretes onto the protoplanet, while the remainder accretes 
onto the circumplanetary disk. The gas accreting onto the circumplanetary disk orbits in 
the disk and gradually falls into the protoplanet. The density contours in wall surfaces show 
a thick torus-like configuration of the disk in model M005A2, and a very thin disk in model 
M1A2. 

Figure [9] shows the mass flux (p|fr|) for model M005A2, in which only mesh points 
having f ^ < are plotted in color. Thus, the gas in the black region has a positive radial 
velocity Vr > and flows out from the protoplanetary system. Figures ^jp, and b indicate 
that a large part of the gas enters into the protoplanetary system from the vertical direction 
around the Hill sphere. However, the gas does not enter into the protoplanetary system on 
the equatorial plane. Figure |9]c shows that, even when the gas has a negative radial velocity 
Vr < inside the Hill sphere (f < fn), it cannot enter into the region r < r}i/2 on the 
equatorial plane: the flow turns round and returns (see, for example. Figure [7]). Therefore, 
althou gh the gas enters the Hill s phere , it cannot reach the protoplanet on the equatorial 



plane. iPaardekooper fc Mellema] ( 120081 ) also showed that the gas flows out from the Hill 



sphere on the equatorial plane (the equatorial outflow). 



Machida et al.l (|2008[ ) showed that the specific angular momentum of the protoplanetary 
system increases in proportion to j oc Mp. They also showed that the circumplanetary 
disk appears when the gravitational energy exceeds the thermal energy in the whole region 
of the Hill sphere, which occurs when Mp > 0.08Mjup at the Jovian orbit. Thus, when 
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the protoplanet mass exceeds Mp > O.OSMjup at the Jovian orbit, the circumplanetary 
disk appears, and then it increases its size and mass with the protoplanet mass because 
the massive protoplanet acquires a larger amount of angular momentum from the shearing 
motion in the protoplanetary disk. One can expect that the emergence of the circumplanetary 
disk may reduce the accretion rate. However, we confirmed that the circumplanetary disk 
keeps a constant mass, indicating that the circumplanetary disk does not strongly affect the 
gas accretion rate onto the protoplanet. In summary, although the flow pattern around a 
massive protoplanet is different from that around a less massive protoplanet, fine structures 
such as the circumstellar disk around the protoplanet does not strongly affect the accretion 
rate onto the protoplanet. Thus, once the gas is gravitationally captured by the protoplanet, 
it accretes onto the protoplanet. 



4.5. Thermal Effects on the Mass Accretion Rate 

As shown in FigureEl when the protoplanet is less massive (M < 0.2Mjup), the accretion 
rates for models with different equations of state have almost the same value. On the 
other hand, when the protoplanet mass exceeds Mp > 0.2Mjup, accretion rates are slightly 
different for models with different equations of state. This slight difference is considered to 
be caused by thermal evolution around the protoplanet. The deviation from an isothermal 
approximation is sm all in a le s s mas sive system, while it becomes slight large in a massive 



system, as shown in iMachidal (120091 ). The density distribution around the protoplanet for 
models Mil (left) and MlAl (right) is shown in Figure [TOl in which the circumplanetary 
disks are represented by constant density surfaces. The figure indicates that the isothermal 
model (model Mil) has a thinner disk than the adiabatic model (model MlAl). The orange 
constant density surface of p = 10^ shows that model MlAl has a torus-like disk that flared 
up outwardly, while model MlAl has a thick disk. The red constant density surface of 
p = 10^ shows that, near the protoplanet, model MlAl has a more compact disk than model 
MIL 

Figure [TT] illustrates the mass of the gas envelope (or the circumplanetary disk), which 
is integrated from the centre. The flgure indicates that a model with a harder equation of 
state has a more massive envelope. Since the circumstellar disk in these models is sufficiently 
gravitationally stable (Toomre's Q S> 1), the gas in the circumplanetary disk barely accretes 
onto the protoplanet with dynamic instability. However, as shown in Figure [TUl the model 
with large thermal energy, that is, a harder equation of state, has a more massively flared 
disk. In such a disk, the effect of centrifugal and gravity forces is relatively small, because the 
thermal energy is relatively large. Therefore, it is likely that the accretion rate in such a disk 
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becomes higher than that in a model with lower thermal energy in which the circumstellar 
disk is strongly supported by the centrifugal force, because the path of the gas streamlines 
accreting onto the protoplanet is changed. Finally, it follows that the massive disk has a 
larger accretion rate. 

At last, we comment on a mass fraction of the adiabatic gas. For less massive proto- 
planets, owing to the shallow gravitational potential, the density enhancement is small even 
near the protoplanet. Since the gas density is connected to the gas temperature through the 
equation of state (eq. [T]), such models may have a small fraction of the adiabatic gas. We 
estimated a mass fraction of the adiabatic gas for model with fn = 0.29 (Mp = O.OlMjup) 
and Peri = Po (model MOOlAO). The mass fraction / of the adiabatic gas to the total mass 
within the Hill (r < th), half Hill (r < l/2rn) and 1/10 Hill (r < l/lOrn) radii are / = 0.26, 
0.55 and 1.0, respectively. In addition, the fraction inside the Bondi radius (r < tb = 0.25rH) 
is / = 0.69. As shown in §4.3[ the Bondi radius is more important to investigate the gas 
accretion when re < th. A large fraction of total mass is adiabatic gas inside the Bondi (or 
Hill) radius, while there is small difference in the mass accretion rate among isothermal and 
adiabatic models, as seen in Figure [5l As a result, we concluded that the mass accretion 
rate barely depends on the thermodynamics around the protoplanet. 



5. Comparison with Previous Simulations 



was investigatec 


in some previous studies ( 


Kiev et al. 


2001; 


/. 1 . 1 . 

Tanieawa & Watanabe 


2002 




D'Aneelo et al. 


2003 


Bate et al.ll2003 


ilKlahr & Klevlbood 


Dobbs-Dixon et a 




2007 




I'aardekooDer & Me 


2008; 


Fouchet & Maver 


20081: 


D'Aneelo & Lubow 


20081; 


Avhffe & Bate 


2009 


) . Subsections 



§5. II to [5^ will explain the classification of the previous simulations by classifying them into 
three categories and discussing the salient features of each category. 



5.1. Two-dimensional Isothermal Calculations 



Tanigawa fc Watanabd ( 12002| ) derived the mass accretion rate under the isothermal gas 



condition in a two-dimensional local simulation. The spatial resolution of their calculation 
is comparable to our study. Howe ver, the mass accretion r ate is greatly different from the 
results obtained in this study. In iTanigawa fc Watanabd ( 12002| ). the mass accretion rate 
continues to increase as a function of the protoplanet mass, while it saturates at Mp ~ 
0.2Mjup in our three-dimensional calculation. Furthermore, although in this study, the 
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mass accretion rate increases wit h the protoplanet mass i n the range of Mp < 0.2Mjup, 
its value is much smaller tha n in iTanigawa fc Wataiiabd (120021) . For example, at Mp = 
O.lMjup, the accretio n rate in ITanigawa fc Watanabd (I2OO2I ) is M = 2.3 x 10"^Mjupyr~^ 
(see, Equation [19] of Ixanigawa fc Watanab3 l2002l ). while it is M = 2.7 x lO'^Mjupyr-^ 
in our isothermal model. Thus, in the range of Mp < 0.2Mjup, the accretion rate based on 
a two-dimensional calculation is two orders of magnitude larger than that obtained from a 
three-dimensional calculation. The difference becomes larger in the range of Mp > 0.2Mjup, 
because the accretion rate continues to increase in two- dimensional c alcula tions, while it 
gradually decreases in the three-dimensional calculations. iKley et al. I fl200ll ) compared the 
mass accret ion rate in their three-dimensional simulations with that in two-dimensional 
simulations (jKleylll999l ) and showed that the mass accretion rate d erived in three-dimensional 
simula tions is smaller than that in two-dimensional simulations. iPaardekooper fc Mellema 
(120081 ) also commented on the decrease in mass accretion observed in three-dimensional 
simulations compared with two-dimensional simulations. The difference in mass accretion is 
considered to be caused by the spatial dimensions. 



5.2. Three-dimensional Isothermal Calculations 



Kley et al. 



(120011 ). iD'Angelo et al. I (120031 ). and iBate et al.l (l2003f ) estimated the mass 
accretion rate under the locally isothermal assumption in global simulations. The mass ac- 
cretion rates derived from these global simulations corresponds well to the rate derived from 
our local simulations, despite the fa ct that the spatial reso lution and size of the sink are 
considerably different. For example, iD'Angelo et al. I (120031 ) estimated the mass accretion 
rate in their global three-dimensional simulations with a spatial resolution of Ax = 0.06 fn, 
which is 40 times coarser than ours, and defined the sink as the region inside one-tenth 
of the Hill radius {rsmk < 0.1 fn), which is 28 times large r than our sink s i ze (se e §3.3p . 



Thus, our study can resolve the present Jovian radius, while iD'Angelo et al. I (120031 ) cannot 
resolve the radius. As well, in their calculation, the circumplanetary disk cannot be suf- 
ficiently resolve d because the circumplanetary disk i s formed in the region r < 10 — 50 rp 
(lMachidall2009l ). In summary, the differences between iD'Angelo et al. I (120031 ) and our study 
are the spatial resolution and numerical setting. We calculated the evolution of the proto- 
planet aryjy;stem_jn__th^lo^^ with finer spatial resolution, while this cannot be 
done in iD'Angelo et al. I (j2003l )'s global simulation with coarser spatial resolution. Despite 
the differences, the accretion rate in our isothermal models corresponds well to their results 
(compare Figure [5] with Figure 7 of iD'Angelo et al. II2003I ). The mass accretion rate has a 



peak value of M ~ 8 x lO'^Mjupyr"^ at Mp = 0.2 
while it has a peak value of M ~ 4 x 10~^Mjup yr^^ 



0.3Mjup in iD'Angelo et al. I (|2003i), 
0.2Mjup in our study. Fur- 



at Mp ~ 
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thermore, the mass accreti on rates at M = 0.01 and IMjup are M ~ 3 x 10 '''Mjupyr ^ 
and ~ 4 X 10~^Mjup yr"-"^ in iD'Angelo et al. I (120031), while they are 8.2 x 10~^M.i„p yr"^ and 



3.6 X 10 ^Mjupyr ^ in our study. It should be noted that although iD'Angelo et al. I (120031 ) 



investigated the evolution of the protoplanetary disk by changing the gravitational potential 
of the protoplanet, the mass accretion rates differed slightly in each model. 



Kley et al. I (120011 ) investigated the evolution of the protoplanetary disk only when the 



protoplanet has a Jovian mass (IMj^p) under the isothermal gas conditio n in their three- 



dimen sional simulations, in which the spatial resolution is coarser than iD'Angelo et al. 



(120031 ). They found that a mass accretion rate of 6 x 10~^Mjup yr~^, while it is (2.2 — 4.1) x 
10~^yr~^ at IMj^p in our study ( see, f or example. Fig. [5]). Thus, our results correspond 
well with their results. iBate et al.l (l2003l ) also calculated the evolution of the protoplanetary 
di sk in a three-dirnension al global simulation with a sp atial resolut i on co mparable to that 
of iD'Angelo et al. I (120031 ). The accretion rate derived in lBate et al.l (120031 ) corresponds well 
to ours within a factor of three. In addition, it shows the same trends as in our study. That 
is, the accretion rate has a peak at Mp ^ O.lMjup and decreases for Mp > O.lMjup. 

In summary, all three-dimensional calculations show almost the same mass accretion 
rate, though they have adopt different calculation settings (i.e., local or global calculation), 
spatial resolutions and sizes of the sink radius. In these calculations, the Hill radius or the 
region near the Hill sphere is resolved, while the region near the protoplanet, much smaller 
scale than the Hill radius is not always resolved sufficiently. This indicates that the mass 
accretion rate is regulated by the region around the Hill sphere, not by the details at smaller 
scale length. Thus, we conclude that the mass accretion rate is safely estimated when the 
Hill radius is resolved with adequate spatial resolution. 



5.3. Three-dimensional Radiative Calculations 



Recentlv.lKlahr fc KlevI (120061 ) . IPaardekooper fc Mellemal (120081 ) . iFouchet fc Maveii (120081 ) 
and lAyliffe fc Bate! (120091 ) investigated the evolution o f the protoplaii e tary d isk using three- 
dimensional, radiation-hydrodynamical simulations. iKlahr &: Kleyl (120061 ) estimated the 
mass accretion rate of 5.1 x 10~^Mjupyr~^ when the protoplanet has a Jovian mass. They 
concluded t hat the mass accretion rate in a r adiative model is larger than that in an isother- 
mal model. iPaardekooper fc Mellemal (120081 ) calculated the evolution of the protoplanetary 
disk with several models with different protoplanetary mass under the locally isothermal 
approximation, while they calculated it with the radiation-hydrodynamical simulations only 
with protoplanet mass of 0.6 M® and The accretion rates in their isothermal calcula- 

tion are identical to our results. As well, the accretion rate in their radiation-hydrodynamical 
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simulations is also comparable to our results. It should be noted that, since the range of 
the protoplanetary mass in their radiation-hydrodynamical simulations is different from our 
study, a direct comparison of the results cannot be performed. They concluded that al- 
though the accretion rate in the radiation-hydrodynamical simulation is smaller than that 
in the isothermal simulation, the difference is not dramatic. It should be noted that, for 



radiation-hydrodynamical simulations, iPaardekooper fc Mellemal (120081 ) only calculated the 
disk evolution with very small protoplan etary masses of 0.6 M m and SMq. These masses are 
considerably different from that used by iKlahr fc Kleyl (120061 ). 



Fouchet fc Mayerl (12008! ) investigated the evolution of the protoplanet disk including 
both the self-gravity and radiation physics for a protoplanet with a Jovian mass, IMjup. They 
also calculated the protoplanetary disk under the isothermal approximation and compared 
them with the radiation-hydrodynamical simulation. They obtained an accretion rate of 
5 X 10~^Mjupyr~^ in the isothermal simulation and ~ 2 x 10~^M.T„n vr~^ in the radiation- 



'.Tnp yr 



hydrodynamical simulation. These rates are well identical to iKlahr fc Kiev! (120061 ) and our 



Avliffe fc Batel J2009h also estimated the 



results of (2.2 - 4.1) x 10"^ yr'^ at Mp = IMjup. 
mass accretion rate in their radiation-hydrodynamical simulation and found an accretion 
rate similar to those derived in previous studies. 



Howe ver, the radiative ef f ect on the mass accretion rate is controversial. iPaardekooper fc Mellema 



(12008! ) and lFouchet fc Mayerl (120081 ) showed that t he accretion rate in a noniso thermal disk is 



reduce d compared to the isothermal disk (see also. lAyliffe fc Batd!2009l ). while iKlahr fc Kley 



(I2OO6! ) showed the mass accretion rate in radiative model is larger than in the isothermal 
model. In our study, the accretion rate for the isothermal model is smaller than for the 
adiabatic models. However, the difference in the accretion rate between the isothermal and 
radiative models is not so large. This can b e attributed to the gravitational energy of the 



protoplanet dominating the thermal energy (!Machidal!2009l ). Thus, it is expected that the 



accretion rate barely depends on the thermal evolution in the protoplanetary disk. For ex- 



range (2-6) x lO-^Mju^yr-^ ( 
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8). This indicates that 



the growth timescale of the Jovian planet at Jovian orbit is ~ 10^ years. 



5.4. Disk Viscosity and Effect of the Gap on the Mass Accretion Rate 



In this subsection, we discuss the relation between gap formation and the mass accretion 
rate onto the protoplanet. Some previous studies pointed out that the redu ction of the mass 
accretion rate for massive protoplanets is related to the gap formation (e.g.. iD'Angelo et al. 
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20031 : iBate et al.l l2003l ). The density gap in the protoplanetary disk is the result of the 
competition between torques exerted on the disk by the planet and by the disk itself. The 
planet gives angular momentum to the outer part of the disk, and it takes angular momentum 
from the inner part of the disk (iGoldreich fc TremaindllQSOl ). Thus, a protoplanet tends to 
open a gap. On the other hand, the disk kinematic viscosity makes the gas back to the gap. 
As a result, t he gap formation depends on both the disk kinematic viscosity and mass of the 
protoplanet ( iLin &: Papaloizoulll986l : iBryden et al.lll999l ). A massive protoplanet in the disk 
with a smaller viscosity tends to show a clear gap. However, the formation condition, size, 
and width of the gap are not clearly understood. 

In the inviscid (z/ = 0) gas disk model as adopted in this study, it is considered that 
the gap formation o ccurs when the Hill radius exceeds the scale height of the protop lanetary 
disk, i.e., > h ( jLin fc Papaloizoul Il986l : iBryden et al.l Il999c ICrida et al.l 120061 ). Using 
equations and (|9]), this condition can be described as 



M 



J up 



^ 0.12 (- 
^ VI 



AU 



3/4 



(17) 



which indicates that, at Jovian orbit, the density gap begins to appear when the mass of 
the protoplanet exceeds Mp > 0.4Mjup. As shown in §4.31 we found that the gas flows 
into the planetary system in the region x = 1 — 3 th when the protoplanet is more massive 
than Mp > 0.08 Mjup at Jovian orbit. Thus, when the gap width Agap becomes larger 
than Agap ^ 2 — 6rH, the mass accretion rate onto the protoplanet is expected to decline. 
The decline of the mass accretion rate at Mp > 0.3 — IMjup in global simulations (e.g.. 



D'Angelo et al. 1 12003|) is caused by a wide gap formation. On the other hand, our local 
simulation shows no clear gap, though the low-density region is shallower than that in the 
global simulations. Local simulat ions are not appropr i ate to treat gap formation b ecause 
of the radial boundary condition ( Miyoshi et al. Ill999l : iTanigawa &: Watanabell2002l ). since 
the gap properties, such as the gap depth and width, depend on the size of the simulation 
box. Therefore, when the mass of protoplanets exceeds Mp ^ 0.4Mjup at Jovian orbit, it is 
expected that the mass accretion rate derived in the local simulation is smaller than that in 
the global simulation owing to the gap. On the other hand, in the range of Mp < 0.4Mjup, 
the mass accretion rate derived in the local simulation is always applicable, because no clear 
gap forms in this mass range. In addition, the mass accretion rate derived in this study 
would be valid even for Mp > 0.4Mjup when the disk viscosity is sufficiently large, because 
no clear gap appears in such a viscous disk. Moreover, the mass accretion rate in the local 
simulation well agrees with that in the global simulation in the range of Mp < IMjup as 
denoted in §5.21 This may be because the density gap barely affects the gas accretion in this 
mass range. However, global simulations show a rapid decline of the mass accretion rate for 
Mp > IMjup, in which the density gap seems to significantly affect the mass accretion. In 
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such a case, the mass accretion rate derived in the local simulation cannot be applicable. 

In this study, the mass accretion rate is derived in the local simulation. This accretion 
rate may be overestimated when the protoplanet is sufficiently massive and the protoplan- 
etary disk has a sufficiently small viscosity, because the gap formation decreases the mass 
accretion rate. Thus, the mass accretion rate obtained here corresponds to the possible 
maximum value that is attainable without forming a density gap in the disk. However, even 
when the density gap is formed, the formula of the mass accretion rate is applicable with 
the reduced disk surface density owing to the gap, because we investigated the gas flow in a 
dimensionless form, and all physical quantities are scalable, as denoted in Equation f|T3l) . 



6. Growth Timescale of Gas Giant Planets 



As mentioned in ^ to form a gas giant planet, the solid core with a mass of 1 — 
10 Mf^ needs to acquire th e gas component from the protoplanetary disk. At the end of the 
"oligarchic growth" stage (IKokubo fc Idalll998l ). several protoplanets composed of the solid 
component, which will be referred to as the planetary core, appear with an orbital separation 
of about 10 Hill radii f lKokubo fc Idall2002[ ). Only when a planetary core acquires sufficient 
gas from the protoplanetary disk before the dissipation of the protoplanetary gas disk, the 
gas giant planet is possible to form. The mass and growth timescale of the planetary core 
were investigated in detail using A^-body simulations. In the remainder of this section, using 
the mass accretion rate onto the gas giant planet obtained in this study, as well as the mass 
a nd g rowth timescale of the planet core obtained from A^-body simulations (IKokubo &: Ida 



2OO2I ). an investigation of the final mass and growth timescale of the gas giant planet is 
performed. 



Kokubo fc Idal (120021 ) derived an (isolation) mass and growth timescale of the planetary 
core as a function of disk model. Adopting b = 10 (the orbital separation normalized by 
the mutual Hill radius of the cores), a = 3/2 (the power index of the radial surface density 
distribution), and = Mq (the mass of the central star) in Equations (17) and (26) of 
Kokubo &: Idal (120021 ). the isolation mass Miso and growth timescale Tcorc of the planetary 



core can be described as 
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where /ice (ice factor) and Si (reference solid surface density at 1 AU) can be related to the 
solid surface mass density as 



^soiid = /iceSi[^j gcm-3. (20) 



The ice factor /ice expresses the increase of solid by ice condensation over the snow-line agnow 
and /ice = 1 (a < asnow=2.7 AU) and 4.2 [ a > asnnw)- As the re ference solid surface density, 



we adopt the minimum-mass disk model flHayashi et al. Ill985l ) with Ei = = 7. In the 



following, the solid surface density is parameterized as Si/S^ = 1/2, 1 and 2. We set the 
mass of the planetary core as the isolation mass Mcore = 



'■ISO- 



In a conventional scenario of a gas planet formation (e.g.. lMizundll980l : IStevensorull982 



Bodenheimer fc Pollackl Il986l ). the planetary core is heated by the continuous accretion of 
planetesimals, and has a hydrostatic gas envelope when a planetary core is less massive 
than Mcore ^ 10 M®. Then, the planetary core cannot sustain the gas envelope and the 
rapid gas accretion is triggered (i.e., the runaway gas accretion phase) aft er it grows up 



to Mcnre x^, 10 Mq. However, the current planetary accretion theory (e.g., iKokubo fc Ida 



19981 . |2002[ ) indicates that the accretion of planetesimals almost halts after the isolation core 
formation (or after the oligarchic growth stage). In this case, the planetary core cannot have 
a hydrostatic envelope owing to the absence of the heat source even when the planetary core 
is Mcore ^ 10 M®. Therefore, the runaway gas accretion is triggered just after the isolation 
core formation. 



Figure [T2]a shows the mass of the planetary core for different Si based on Equation ( fTSl) . 
In the figure, the sudden rise of the core mass at 2.7 AU is caused by the increase of solid 
component at the snow-line. In the standard model (Si = S^, the red line), the planetary 
core has a mass of 2.8 M® at the Jovian orbit (5.2 AU) and 4.4 M® at the Saturnian orbit 
(9.6 AU). 

Figure [T2b shows the growth timescales of the planetary core Tcore (thin solid lines. 
Equation [19]) and Jovian-mass gas planet Tgas (thick solid lines) in each orbit. Using Equa- 
tion (fT2|) . the mass accretion rate is integrated as a function of time until the planetary core 
reaches the Jovian mass Mjup in each orbit using 

Mj^p = M,,,, + J -^dt, Mp(t = 0) = Mcore, (21) 

where L = Lq, T = Tq and po = Po(Si/S^) are used in Equation ( fT2l) . Thus, Tgas represents 
the time required to form a Jovian-mass gas planet after the formation of the isolation core. 
The figure indicates that Tg^s ^ ^core in any orbit. For example, in the standard model, 
that is. Si = S^, the growth timescale of the planetary core is Tcorc = 2.7 x 10'' yr at the 
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Jovian orbit and 1.7 x 10^ yr at the Saturnian orbit, while that of a Jovian gas planet is 
4.7 X 10^ yr at the Jovian orbit and 7.7 x 10^ yr at the Saturnian orbit. Therefore, 



'gas 



the gas-accretion timescale is two or three orders of magnitude shorter than the core growth 
timescale. 



As shown in lKokubo fc Idal (120021 ). gas giant planets can form in the limited range of the 
protoplanetary disk where the planetary core is large enough to capture the gas component 
before the dissipation of the gas disk. The solid curves of Figures [T2]c and d show the planet 
mass for the disk lifetime of tgas = 5 x 10^ yr (c) and 10^ yr (d). In these figures, using 
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Equation f[T2|) . the planet mass Mp at each orbit is derived by 



dM, 



dt 



^dt, Mp(t = 0) = Mc, 



The dotted lines of Figures [T2b a nd d indicate the max imum mass in the ring of the proto- 
planetary disk around the planet (IHayashi et al. Ill985l ). which is described as 



27rap 6r Eg 



(23) 



where 6r is the ring width and is expressed by 6r = 2/HrH, where /h is the parameter 
representing the ring width. Adopting Mp = Mmax and M^. = Mq in equation ([9]), 6r can 
be described as Sr = 2/h ap(Mmax/3 Mq)^/^. Thus, Equations (123!) gives 



M„ 



-1/2 



(24) 



In equation (124|) . Eg = Sg is adopted as the gas surface density of protoplanetary disk, 
where /g is related to the solid surface density as /g = Si/S^ because the constant ratio of 
the gas to solid surface d ensity is assumed. Th e fiducial gas surface density Eg is given by 
the standard disk model (IHayashi et al. Ill985l ) as 



= 1.7 X 10^ 



1 AU/ 



-3/2 



gem 



(25) 



Equation is the mass of the ring in the protoplanetary disk in the range of < fnrn, 
where the origin of x is the position of the planet. As shown in §4.31 our calculation indicates 
that the gas in the range of |x| ^ 3rH flows into the protoplanet system. Thus, /h = 3 is 
adopted in equation (l24ll . Since the protoplanet cannot acquire the gas component exceeding 
the mass in this ring, the mass M^ax means the maxir num mass of the ga s component 
acquired by the protoplanet at each orbit (for details see. IHayashi et al. 1119851 1. 



The thick gray curves of Figures lT2b and d are the planet mass attainable at each orbit. 
Figure lT2lc indicates that the planet mass is determined by the dotted line (i.e., eq. [23]) for 
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flp < 5AU, while by the sohd curve (i.e., eq. [22]) for Op > 5AU. For < 5AU, since the 
growth timescale of the gas planet is shorter than the lifetime of the gas disk, the planetary 
core accumulates all gas in the ring with < Srn in t < 5 x 10^ yr. On the other hand, 
for Op > 5AU, since the gas disk dissipates before the protoplanet accumulates all gas in 
the ring, the planet mass is determined by equation (1221) . As shown in equation (fT2l) . since 
the mass accretion rate is a decreasing function of the orbital radius ap, the planet mass 
decreases with Op for Op > 5AU. The thick gray curve in Figure [T2b shows that the planet, 
respectively, has almost the Jovian and Saturnian mass at Jovian and Saturnian orbit. For 
ttp > 10 AU, the mass of the gas component is less massive and comparable to the solid core 
mass (Mp ~ Mcorc ~ 10 M^). In addition, t he gas planet forme d at a„ > 2.7 AU has a core 
mass of IM® < Mcore < 10 M®. Note that ISaumon &: GuillotI (120041 ) studied the internal 
structure of Jupiter and Saturn, and estimated the core mass of the Jupiter in the range of 
— 11 Mq, and core mass of the Saturn in the range of 9 — 22 M^. These features (mass 
of gas and solid components at each orbit) well correspond to those of giant planets in our 
solar system (Jupiter, Saturn, Uranus, and Neptune). 

Although we assumed that the gas disk dissipates in t = 5 x 10"^ yr after the isola- 
tion core formation in Figure \T2b. observations indicate that the protoplanetary disk has 
a hfetime of ^ 0.3 - 30 M yr Jflaisch et all I2OO1I : iHartmannl I2OO5I : Isilverstone et all I2OO6 



Flaherty Sz MuzeroUd l2008l ). Figure \T2H shows the planet mass when the gas disk exists 
for t = 10^ yr after the isolation core formation. In this figure, the planet mass is deter- 
mined only by the dotted line (i.e., the maximum mass in the ring, eq. [M])- As shown in 
Figure [T^. since the growth timescale of a gas planet is shorter than 10^ yr at any orbit, 
the protoplanet can accumulate all gas component around it when the gas disk exists for 
t > 10^ yr. The thick gray line in Figure [T2] g? shows that although the Jovian mass planet 
is formed at Jovian orbit, more massive planets than Jovian mass are formed at Saturnian, 
Uranian, and Neptunian orbits, which contradicts our solar system. 

Figures [T2b and d indicate that the gas giant planet with the core mass of Mcorc > 1 
can be formed in the orbital range of Op > 2.7 AU. They show that the mass of planets 
exceeds Jovian mass in the range of Op > 5AU when the gas disk exists for tg^s > 10^ yr 
after the isolation core formation. Thus, to realize our solar system in the present model, it 
is expected that the gas component of the protoplanetary disk needs to dissipate just after 
the isolation core formation. 

We need to pay attention to viscous evolution of the disk for estimating the planet 
mass. In Figures [T^lc and d, the planet mass in the range of Op < 2.7 AU for disk with 
Tgas = 5 X lO'^yr and in whole orbital range for disk with Tgas = 10^ yr is determined by the 
disk surface density at the location of the planetary core formation when the gas accretion 
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timescale is shorter than the disk hfetime. In other words, the planet mass is determined 
by in situ surface density (or in situ disk mass). However, in these figures, the disk viscous 
evolution is ignored. If the gap is refilled by the viscous evolution, the planet can acquire 
the gas component exceeding Mmax- Thus, the planet mass becomes more massive than gray 
line in these figure, when the disk has a significantly large viscosity. 

At last, we comment on the hydrostatic envelope. As shown in the above, to estimate the 
formation time of a gas giant planet, we assumed no planetesimal accretion onto the planetary 
core after the isolation core formation. Although this assumption is supported by recent 
numerical simulations, the formation timescale of a gas giant planet may be considerably 
different when the planetesimal accretion does not completely halt even after the isolation 
core formation. With the planetesimal accretion, the planet ary core can sustain a hydrostatic 



envelope for ~ 10*^ yr at the maximum ( llkoma et al. 



20001 ). In this situation, the formation 



time (or gas accretion time) for the gas giant planet is prolonged. Thus, Tgas in Figure [12 
gives the minimum time necessary for the formation of gas giant planets. 



7. Summary and Discussion 

In this study, a calculation of gas accretion onto a planetary core including the thermal 
effect was performed, in which the mass accretion rate and growth timescale of a gas giant 
planet were estimated. 1, 

The mass accretion rate derived in this study quantitatively agrees with those derived in 
previous simulations. It is surprising that the same value for the accretion rate is derived from 
three-dimensional simulations with different parameter values, including parameters such as 
spatial resolution, the size of the sink, and the treatment of the thermal effect, including 
locally isothermal approximation, barotropic equation of state, and radiative hydrodynamics. 
This suggests that the mass accretion rate is only determined by 'the protoplanet mass' and 
'the properties of the protoplanetary disk.' In our calc ulation, we r e solve d the Jovian radius 



of the cell width. On the other hand, for example, iKley et al. I (l200ll ) derived the mass 
accretion rate by using the sink radius of rH/2. In §4.41 we showed that almost all gas 
falling into the half Hill radius (r < rH/2) can accrete onto the surface of the protoplanet. 
Although the radiative effect can affect the mass accretion rate, the difference is less than a 
factor of ~ 3. As a result, this study, like previous studies, shows that the mass accretion 
rate can be accurately estimated for resolution higher than half Hill radius. In other words, 
the high spatial resolution that resolve Jovian radius in not necessary to investigate the mass 
accretion rate and growth times scale of gas giant plants. 
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A comparison of our results with the A^-body simulations for the solid core aggregation 
shows that the gas giant planet can form in the orbital range of Op > 2.7 AU in a short 
timescale of ~ 10^ — 10^ yr after the planetary core inside the gas giant planet is formed in 
a longer timescale of ~ 10^ yr. In each orbit, the growth timescale of the gas giant planet 
is two orders of magnitude shorter than the core aggregation timescale. This may indicate 
that, in our solar system, the planetary core with 1 — 10 began to capture the gas 
component to form gas giant planets such as Jupiter and Saturn just before the dissipation 
of the protoplanetary disk was completed. Otherwise, the protoplanet accumulates all gas 
around it before the gas dissipation when the gas disk exists for ^ 10^ yr. In the latter case, 
the mass of the gas giant planet may be determined only by disk properties (i.e., the surface 
density and density profile of the protoplanetary disk). Thus, to investigate the formation 
of the gas giant planet, we need to precisely estimate the properties of the protoplanetary 
disk, and the gas dissipation process (or timescale). 

In this study, it was not possible to investigate gap formation in detail, since the sim- 
ulations were limited in a local region around the protoplanet. However, it was possible to 
resolve the present Jovian radius and the circumplanetary disk. To understand gap forma- 
tion and the mass accretion rate, it is necessary to simulate a protoplanetary system using 
a global simulation that can resolve the protoplanet. However, such a calculation requires a 
large amount of CPU time. As well, issues concerning disk viscosity would need to be better 
understood before further investigation of both the 'protoplanetary' and the 'circumplane- 
tary disk' can be performed. In general, it is considered that disk viscosity is caused by the 
magneto- rotational-instability (MRI). Thus, to properly calculate the viscous disk (circum- 
planetary and circumstellar disks), it is necessary to include magnetic effects. Although it 
is currently possible to calculate t he evolution of the magnetized disk with a lower spatial 



resolution (IMachida et al. Il2006bl ). a simulation with a higher spatial resolution (or more 
CPU power) is necessary to investigate the disk viscosity correctly. For such a simulation, 
the next generation of supercomputers is required. However, the mass accretion rate derived 
in the present local simulation is consistent with those obtained in global simulations. Thus, 
we could link the mass accretion rate in a global simulation to that in a local simulation. 
Therefore, we can safely estimate the growth of the gaseous protoplanet. 
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Table 1: Model parameters 



Model 



m re Mp (5.2 AU) " 



rsi„k(10-3) 



(MOOIAO, MOOlAl, M001A2, M001A3, MOOII) 0.29 0.073 0.01 

(M005A0, M005A1, M005A2, M005A3, M005I) 0.5 0.38 0.05 

(MOIAO, MOlAl, M01A2, M01A3, MOII) 0.63 0.75 0.1 

(M02A0, M02A1, M02A2, M02A3, M02I) 0.8 1.54 0.2 

(M04A0, MOlAl, M04A2, M04A3, M04I) 1.0 3.0 0.4 

(M06A0, M06A1, M06A2, M06A3, M06I) 1.15 4.56 0.6 

(M08A0, M08A1, M08A2, M08A3, M08I) 1.26 6.0 0.8 

(MIAO, MlAl, M1A2, M1A3, Mil) 1.36 7.55 1 

M02ISS 0.8 1.54 0.2 

M02ISM 0.8 1.54 0.2 

M02ISL 0.8 1.54 0.2 
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Fig. 1. — The density-temperature diagrams for the isothermal (pcri = oo) and adiabatic 
(Pcri =P0) lOpo, 10^ po, and 10'^ po) models. The right and upper axes indicate the dimensional 
temperature and density at the Jovian orbit. 
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Fig. 2. — The mass accretion rate as a function of the sink radius for models M02ISS, 
M02ISM, M02ISL, and M02I. 
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Fig. 3. — The density distribution {color scale) on the cross section in the z = plane 
for models M001A2, M005A2, M01A2, M04A2, M08A2, and M1A2. The dashed circle 
represents the Hill radius. Four grid levels are shown in each panel (/ = 2, 3, 4, and 5). The 
level of the outermost grid is denoted in the top left corner of each panel. The elapsed time 
t is denoted above each panel. 
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Fig. 4. — The gas structure around the protoplanet in a bird's-eye view for models M001A2, 
M005A2, M01A2, M01A2, M04A2, M08A2, and M1A2. Density distribution on x = 0, ?/ = 0, 
and 2 = plane are projected onto each wall surface. The color surfaces indicate constant 
density surfaces: p = 10^ po (red) and 10^ po (orange). Velocity vectors (arrows) are plotted 
on the bottom wall. The size of the domains is shown in each panel. The grid level is shown 
in the top left corner of each panel. 
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Fig. 5. — The mass accretion rate for models with pcri = Po(x)) lOpo (^)) 10^ po (O)) lO^po 
(□), and oo (+) as a function of the cubed Hill radius. The right and upper axes indicate 
the dimensional mass accretion rate and the protoplanet mass at the Jovian orbit. The solid 
line represents the fitting formula. 
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Fig. 6. — The gas streamlines and mass flux p\v\ for models M001A2, M005A2, M01A2, 
M04A2, M08A2, and M1A2. Each top panel: The streamlines are plotted in three- 
dimensions. The color of the streamlines indicates the vertical component of the velocity 
at each mesh point. The values of the Bondi, doublewide Bondi, Hill and doublewide Hill 
sphere are plotted. The grid level is shown in the top left corner of each panel. Each bottom 
panel: The mass flux in the x-z plane is plotted on the bottom ^/-boundary. The gas flows 
into the protoplanetary system are given only in the gray region. 
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Fig. 7. — The gas streamlines integrated from the region inside the sohd square for model 
M04A2. The color of streamlines represents the mass flux at each mesh point. Density 
distribution on the x = 0, y = 0, and z = planes are plotted on each wall surface. The 
distribution of p\v\ on y = —6 h plane is also plotted. The color surfaces indicate a constant 
density region of p = 0.5. 
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Fig. 8. — The gas structure around the protoplanet for models M005A2 and M1A2. Density 
distribution on the x = 0, y = 0, and z = planes are projected on each wall surface. 
The color surfaces indicate high-density surfaces. The streamlines {red lines) are inversely 
integrated from the proximity of the protoplanet. The size of the domains is shown in each 
panel. The grid level is shown in the top left corner of each panel. 
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Fig. 9. — The mass flux (pl^rl) around the protoplanet on the y = (top left), x = 
(bottom left) and z = (right) planes for model M005A2. The gas in the colored region has 
a negative mass flux p\vr\ < 0, and flows into the protoplanetary system, while the gas in 
the black region has a positive mass flux p\vr\ > and flows out from the protoplanetary 
system. The dotted circle indicates the Hill radius. 
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Fig. 10. — Same as Figure IH but for model Mil (left) and MlAl (right). 
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Fig. 11. — Cumulative mass integrated from the center for models MlAl, M1A2, M1A3, 
and MIL The right and upper axes indicate the dimensional mass and radius normalized by 
the Jovian orbit. 
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Fig. 12. — The isolation core mass, or mass of the planetary core, Mcore, (b) growth timescale 
for core Tcorc and Jovian-mass planet Tgas, (c) Planet mass Mp when t^eh = 5 x 10^ yr and (d) 
when Tneb = 10^ yr against the orbital radius Cp. The dotted lines in panels c and d indicate 
the maximum gas mass acquired by the planet at each orbit. The open diamonds and square 
in panels c and d indicate the orbital radius when the core mass exceeds M > IM^ (filled 
diamond) and M > lOM^ (filled square). Each open circle indicates the mass and orbital 
radius of Jupiter and Saturn. 



